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Wave-vector resolved radio frequency (rf) spectroscopy data for an ultracold trapped Fermi gas 
are reported for several couplings at T c , and extensively analyzed in terms of a pairing- fluctuation 
theory. We map the evolution of a strongly interacting Fermi gas from the pseudogap phase into 
a fully gapped molecular Bose gas as a function of the interaction strength, which is marked by a 
rapid disappearance of a remnant Fermi surface in the single-particle dispersion. We also show that 
our theory of a pseudogap phase is consistent with a recent experimental observation as well as with 
Quantum Monte Carlo data of thermodynamic quantities of a unitary Fermi gas above T c . 
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While the existence of a high-temperature superfluid 
phase in the BCS-BEC crossover of a strongly interacting 
Fermi gas is experimentally well established, important 
questions remain as to the nature of the gas above the 
superfluid transition temperature T c . In particular, the 
question of whether or not a pseudogap state exists and 
how to identify it is of importance [l| . This is a question 
that may have relevance to the controversy surrounding 
the pseudogap state in the high-T c cuprates. While the 
origin of this state in the cuprates is a hotly debated 
topic, with atomic Fermi gases we can answer the sim- 
pler question of whether or not strong interactions and 
pairing fluctuations alone can lead to a pseudogap phase. 
This, in turn, tells us whether using such an approach to 
explain the pseudogap phase in the cuprates is a viable 
option or if other mechanisms are required. 

As a function of increasingly strong attractive interac- 
tions, a Fermi gas exhibits a smooth crossover (called the 
BCS-BEC crossover), from a weakly attractive Fermi gas 
with a superfluid transition explained by conventional 
BCS theory, to a Fermi gas where interparticle attrac- 
tions are so strong that the fermion pairs form molecules 
and the gas is well described as a molecular Bose gas with 
a Bosc-Einstein condensation transition. In the BCS 
limit the phenomena of Cooper pairing and superfluidity 
occur simultaneously at the phase transition, while in the 
BEC limit pairing and Bose condensation are decoupled 
with pairing of fermionic atoms into molecules occurring 
well above the condensation temperature. The pseudo- 
gap phase refers to the normal state of a strongly inter- 
acting Fermi gas in the center of this crossover, where it 
is proposed that pairs exist above the superfluid transi- 
tion in analogy with the normal state of the gas in the 
BEC limit. However, unlike the pairs in the BEC limit, 
the pairs in the pseudogap state have many-body charac- 
ter with the underlying Fermi statistics playing a crucial 
role, in analogy with the Cooper pairs of the BCS limit. 
A key prediction of theories of the pseudogap phase is 
that there should be a smooth evolution from the many- 



body pairs in the center of the crossover to the molecular 
pairs in the BEC limit [l|, Q and accordingly, in order to 
verify the existence of a pseudogap phase, it is critical to 
examine the evolution of the spectral function from the 
center of the crossover to the molecular limit Q. 

Based on two recent experiments, conflicting conclu- 
sions have been reached about the existence of a pseu- 
dogap state in the strongly interacting Fermi gas. On 
the one hand, thermodynamic measurements [J] have 
been interpreted as well described by Fermi liquid theory, 
without the need for a pseudogap state. On the other 
hand, momentum-resolved rf spectroscopy which 
measures the single-particle spectral function, has been 
interpreted as evidence for a pseudogap state above T c . 

In this work, we present a theoretical investigation 
of the pseudogap regime based on the t-matrix pairing- 
fluctuation approach of Ref.Q, addressing both the 
single-particle spectral function and the thermodynam- 
ics of the g function of interaction strength in 
the BCS-BEC crossover. We find that, in the pseudo- 
gap regime, the single-particle dispersion back-bends at a 
wave vector kr, near the Fermi wave vector hp , indicating 
the existence of a remnant Fermi surface in this strongly 
interacting gas and the importance of Fermi statistics 
to the pairing. As interactions are increased towards 
the BEC limit, kr, disappears rapidly when entering the 
regime of molecular pairing. This picture is supported by 
a comparison of our theoretical results, where we include 
the effects of the trapping potential, with new experi- 
mental data using momentum resolved rf spectroscopy 
to probe the gas for different interaction strengths. In 
addition, we show that the theory also reproduces the 
observed linear behavior in the thermodynamics. 

By the experimental technique introduced in Ref.0], 
excitations of the trapped gas produced by an rf pulse 
are analyzed by time-of-flight imaging to determine the 
wave vector of the excited atoms once the trap has been 
switched off. The new data are presented with an im- 
proved signal-to-noise ratio at the critical temperature 
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FIG. 1. Experimental (circles) and theoretical (full lines) 
EDC for the trap at T c , for several couplings and wave vectors. 



T c , which is accurately determined as the temperature 
where the condensate fraction disappears. We concen- 
trate in the coupling range 0.0 < (/cfof) -1 ^ 1.0, be- 
cause the evolution of interest from the pseudogap state 
to the molecular Bose gas occurs on the positive side of 
the resonance. Here, o,f is the scattering length associ- 
ated with the Fano-Feshbach resonance and kp is given 
by h 2 k 2 F /(2m) = E F = ?iu; (3iV) 1/3 , where h is Planck 
constant, m the atom mass, N the total number of atoms, 
and wo the average trap frequency (we set h = 1). 

Ultracold Fermi gases are peculiar systems, in that 
their intcrparticle coupling can be increased to the point 
when a description in terms of a gas of molecular bosons 
holds, for which a real gap exists in the single-particle 
spectra. This molecular (two-body) physics is of no inter- 
est in the context of the pseudogap, in a similar fashion of 
molecular binding in vacuum being distinct from Cooper 
pairing at finite density in the presence of a Fermi sur- 
face (cf. footnote 18 of Ref.Q). The question then arises 
about what fermionic feature distinguishes the pseudo- 
gap from the molecular phase. We shall find that the 
back-bending of the dispersion curves obtained from the 
single-particle spectral function A{k,oS) (with wave vec- 
tor k and frequency u>) occurs at a wave vector k^ which 
remains close to hp over a wide coupling range even when 
approaching the molecular limit. We refer to this special 
wave vector as Ul because it is reminiscent of the Lut- 
tinger theorem Q , according to which in a normal Fermi 
liquid the radius kp of the Fermi sphere is unaffected by 
the interaction. 

Figure [1] compares the experimental and theoretical 
energy distribution curves (EDC) at T c for five different 
couplings in the window of interest (see Ref.@ for de- 
tails) . We emphasize that the experimental data bear on 
an absolute normalization, in that only the integral over 
wave vector and energy of the EDC curves (and not the 
separate spectra) has been normalized to unity Q. For 
this reason, there is no independent normalization in the 
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FIG. 2. (a) Dispersions and (b) widths of the low-energy EDC 
peak. Experimental data (circles) and theoretical calculations 
for the trap (full lines) are shown for the same couplings of 
Fig. Q3 and compared with the contribution from the radial 
shell with the largest particle number (dashed lines). In the 
left panels the free-particle dispersion k 2 / (2m) is also reported 
for comparison (thin full lines). 



various panels at different k. This renders quite strin- 
gent the comparison with the corresponding theoretical 
calculations, which in turn contain no adjustable param- 
eters. Good agreement results from this comparison. In 
particular, the theoretical calculations well reproduce the 
asymmetry of the experimental curves between positive 
and negative energies, in addition to the peak positions, 
widths and heights (note how the latter change by about 
one order of magnitude from small to large k). Note fur- 
ther the excellent agreement between the theoretical and 
experimental negative energy tails, and the gradual flat- 
tening of the EDC curves for increasing coupling due to 
the increase of intrapair correlations. 

In Fig. [5] the dispersion and full width at half maxi- 
mum of the peak at lower energies are reported over a 
dense set of k values for the same couplings of Fig. [TJ 
and compared with our theoretical calculations. Note 
that a characteristic back-bending is revealed from these 
dispersions [llj ■ This kind of back-bending is typical of 
a BCS-like dispersion, and is associated with the pres- 
ence of a pseu dog ap in a strongly interacting Fermi sys- 
tem 0, [Bl [l2-14|. In addition, the large values of the 
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FIG. 3. (a) Coupling dependence of the Luttinger wave vector 
fcz, for a homogeneous system at T c , according to the theory 
of Ref. 0] (full line) [the value at unitarity from the QMC cal- 
culation of Ref.[17j is also reported (star)]. The inset shows 
the temperature dependence of \zl at unitarity (full line), and 
compares it with those obtained from the temperature depen- 
dence of the chemical potential of the non-interacting (dashed 
line) and interacting (dashed-dotted line) systems, (b) Theo- 
retical (full line) and experimental (squares) coupling depen- 
dence of fci for the trap system at T c . 



widths (which are at least of the order of Ep) and their 
asymmetric behavior between k < kp and k > kp are 
associated with strong deviations from the expected be- 
havior of a normal Fermi liquid (which requires instead 
the quasi-particle widths to be vanishingly small at kp 
[Toll), and confirm the fact that single-particle states in 
this region constitute poor quasi-particles. Large values 
of the widths are not surprising in the context of the pseu- 
dogap physics that results from pairing fluctuations Q. 
Large widths were also obtained by the self-consistent t- 
matrix approach of Ref. [lfjj , which however masked the 
occurrence of a pseudogap near kp. 

It is relevant to discuss how trap averaging affects 
the above results, because different radial shells in the 
trap correspond to different locations in the coupling- 
vs-temperature phase diagram of the homogeneous sys- 
tem. A reasonable hypothesis is that the radial shell with 
the largest particle number (whose radius r max is esti- 
mated to be (0.5 - 0.6)R F where R F = ^^/(mu^)] 1 / 2 
is the Thomas- Fermi radius) contributes most to the total 
signal. The dispersions and widths contributed by this 
shell at r max are represented by dashed lines in Fig. [2J 
which show good agreement with the complete calcula- 
tion. This indicates that both the back-bending of the 
dispersion relations and the associated large widths are 
not an artifact of trap averaging. 

Despite these deviations from the behavior of a nor- 
mal Fermi liquid, in the experimental data and theoreti- 
cal calculations there yet appears a feature which is pre- 
served from the physics of a Fermi liquid. That is the 
Luttinger wave vector kp where the back-bending occurs, 
which is plotted at T c vs (kpap) -1 in Fig. [3J for a homo- 
geneous [panel (a)] and trapped [panel (b)] system. 



Figure [UJa) shows that for a homogeneous system kp 
drops rapidly to zero when (kpap)^ 1 ~ 0.75, where the 
pseudogap in A(k, ui) turns into a real gap and the molec- 
ular limit is reached. Accordingly, we identify the bound- 
ary between the pseudogap and molecular phases where 
this drop occurs. Along this evolution into the molec- 
ular regime, the disappearance of the underlying Fermi 
surface about occurs when the molecular size becomes 
smaller than the interparticle spacing. The existence of 
a remnant Fermi surface with an enclosed volume con- 
sistent with Luttinger theorem was already pointed out 
by ARPES experiments for the pseudogap phase of high- 
ly superconductors [lH, but its importance for delim- 
iting thepseudogap region was not appreciated in that 
context [19( because the interparticle interaction could 
not be controlled. The inset of Fig. EHa) shows the 
temperature dependence of kp calculated for a homo- 
geneous system at unitarity (full line). At high temper- 
atures when the pseudogap closes up, we have identi- 
fied kp as the value where the dispersion of the peak 
at lower energy in A(k,ui) crosses the chemical poten- 
tial Q. This does not contradict our argument that at 
low temperatures the presence of a pseudogap requires 
an underlying Fermi surface, since at high temperatures 
the underlying Fermi surface of a Fermi liquid is not re- 
lated to a pseudogap. The plot also shows the tempera- 
ture dependence of k^o = ^/2m^i°(T) (dashed line) and 
kfj, = \j2n%n(T) (dashed-dotted line), where n°(T) and 
/i(T) are the chemical potentials of the non- interacting 
and interacting Fermi systems, in the order, at the tem- 
perature T. Note that kp about coincides with k^o , while 
fc M is not related with kp. 

Figure M^>) shows the coupling dependence of kp at 
T c for the trapped system, for which the theoretical pre- 
dictions can be directly compared with the experimental 
data (the latter are obtained by a BCS-like fit to the dis- 
persions of Fig. EJa), as explained in Ref.Q). The good 
comparison that results between theory and experiment 
confirms our identification of kp as the relevant quantity 
for identifying the remnant Fermi characteristics of the 
system in the pseudogap phase. 

However, the occurrence of a pseudogap for a unitary 
Fermi gas above T c has recently been questioned, follow- 
ing a result reported in Ref. [3| where a linear dependence 
of the equation of state as a function of [kpT/ /i(T)] 2 (kp 
being Boltzmann constant) was fitted by the Fermi-liquid 
equation of state and then interpreted (2(3| as evidence 
that the Fermi-liquid theory with no pseudogap can de- 
scribe a unitary Fermi gas above T c . To compare with 
the data of Ref. Q and resolve this controversy, we have 
used the theoretical approach of Ref. Q , which contains a 
robust pseudogap associated with a non-Fermi-liquid be- 
havior consistent with the data obtained by momentum 
resolved rf spectroscopy, also to calculate the thermo- 
dynamic properties of a homogeneous system above T c . 
Figure EJa) reports the pressure in the grand-canonical 
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FIG. 4. Thermodynamics of a homogeneous Fermi gas at uni- 
tarity. (a) Pressure vs [fcsT/^(T')] 2 : Experimental data from 
Ref. [J] (circles) are compared with QMC data from Refs.[2l| 
(squares) and [12] (triangles), and with the t-matrix (full 
line). In the inset, the variable [ksT/ p(T)] 2 is transformed 
to (T/Tf) 2 according to the t-matrix. (b) Energy vs (T/T F ) 2 
at fixed density: Experimental data from Ref. [231] (circles) are 
compared with QMC data from Refs.fU (squares) and [22l | 
(triangles), and with the t-matrix (full line). The inset shows 
the density of states per spin (in units of mkF / (2ir) 2 ) for sev- 
eral temperatures in units of T c according to the t-matrix, and 
contrasts it with the non-interacting (n.i.) result, (c) Specific 
heat per particle vs T/Tf obtained from the t-matrix (full 
line) , the experimental data of Ref. [23[ (circles) , and the QMC 
data of Ref.[2lJ (squares) - the dotted line is a guide to the 
eye for the QMC data. The behavior of the non-interacting 
Fermi gas (broken line) is reported for reference ||. 



ensemble vs [ksT/ fi(T)] 2 as in Ref.Q, and shows that 
the linear behavior seen in the experimental data and 
QMC calculations also results from our t-matrix ap- 
proach, both above and below the temperature at which 
the pseudogap appears (indicated by the vertical arrow) . 
The inset of Fig.[4|a) shows that this linear behavior can 
be ascribed to the pronounced temperature dependence 
of the chemical potential, because a non-linear behavior 
results when transforming [k B T / fi(T)] 2 to (T/T F ) 2 over 
the relevant range. The same ch ang e of variables can 
be performed in the experimental |23| and QMC pi. 22| 



data, to obtain the total energy in the canonical ensem- 
ble as a function of (T/T F ) 2 reported in Fig. S^b). This 
shows that in the new variable the linear behavior is lost. 

Yet, it remains difficult to appreciate directly from this 
thermodynamic quantity the presence of a pseudogap in 
a unitary Fermi gas above T c even by the t-matrix calcu- 
lation, despite the fact that a pseudogap is clearly present 
in the single-particle density of states obtained by the t- 
matrix as shown in the inset of Fig.^b) where deviations 
from the non-interacting behavior \J {u + pL{T c ))/Ep are 
evident. Accordingly, by suitable numerical differentia- 
tion of the energy data we have obtained in Fig. @|c) 
the specific heat vs T/Tp. A sharp upturn of this ther- 
modynamic quantity, beginning at a temperature T* well 
above T c where the pseudogap sets in, results clearly from 
the t-matrix calculation, and it is also visible from the 
QMC data at the corresponding value of T c . 

The experimental data in Fig. @|c) appear too scat- 
tered to draw definite conclusions about the presence of 
the upturn and thus of a pseudogap above T c . It should 
be mentioned, however, that a similar upturn of the spe- 
cific heat at a temperature T* above T c was measured in 
underdoped high-T c cuprates and interpreted as revealing 
the onset of the pseudogap regime, whereby a "residual 
superconductivity" remains far above T c [24j. 

In conclusion, we have provided clear experimental and 
theoretical evidence for non-Fermi-liquid behavior in the 
normal phase of a strongly interacting Fermi gas, which 
we have qualified in terms of a pseudogap picture. We 
have further shown that this picture, that appears ev- 
ident in the single-particle dynamics, is also consistent 
with the thermodynamic behavior of the system. 

We acknowledge financial support from the NSF and 
from the Italian MIUR under contract PRIN-2007 "Ul- 
tracold Atoms and Novel Quantum Phases" . 
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Supplemental material: "Evolution of the 
Normal State of a Strongly Interacting Fermi 
Gas from a Pseudogap Phase to a Molecular 
Bose Gas" 

We provide details of the theoretical calculations of 
the wave-vector resolved rf signal and a description of 
the experimental procedures. We also add information 
about the theoretical analysis of the experimental data. 

Pairing-fluctuation theory 

The theoretical approach of Ref. [lj is based on 
a diagrammatic t-matrix approximation, whereby the 
fcrmionic single-particle self-energy S(k, a;) includes pair- 
ing fluctuations. We have used that approach here to 
calculate the single-particle spectral function A(k, to) for 
the homogeneous case: 



4(k,w) =-- 

IT 



Im£(k,w) 



[w-Ck-RcS(k,w)] 2 + [ImS(k,w)] 2 

(1) 

where £k = k 2 /(2m) — /i. For given wave vector k, 
the frequency structure of the real and imaginary parts 
of £(k, w) determines the positions and widths of the 
peaks in ^4(k, w), and is thus responsible for the nontriv- 
ial shape of the dispersions of these peaks vs k = |k|. 

Wave-vector resolved rf spectroscopy 

When final-state effects in the rf transition 0, 0] can be 
neglected (like for the case of 40 K used in the experi- 
ment) , the rf signal in the normal phase is given by 0] : 



RF(cD) 



1 



dr 



dk 



A(k,£(k;r)-cD)/(£(k;r)-d;) 



(2) 

Here, u = w r f — uj a is the detuning frequency where co T { 
is the frequency of the rf photon and u) a the atomic hy- 
pcrfme frequency, r the position in the trap, £(k;r) = 
k 2 /(2m) — yu + V(r) a local energy with trapping po- 
tential V(r) = m (uj x x 2 + uj y y 2 + uj z z 2 ) /2, and /(e) = 
' e e/(k B T) _|_ j\ the Fermi function. [The prefactor in 
Eq.@ is chosen to make the total area of the rf signal 
equal unity.] Equation ([2]) is based on a local-density ap- 
proximation where contributions of adjacent shells in the 
trap are separately considered. 

The rf signal can be analyzed into its individual k- 
components to compare with the experimental technique 
of Ref. 0]. The resulting wave-vector resolved rf signal is 
obtained by dropping the k-integration and considering 
one k-component at a time. More precisely, the selection 
is made over the magnitude fc = |k| while k = k/fc is 
integrated over the solid angle, yielding: 



RF(fc,£>) = 



48fc 2 



drr 2 A(k,£(k;r) 



■a5)/(e(A;;r)-w) 
(3) 
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where the factor k 2 is from the spherical integration and 
r = |r| is the radial position in the trap. The prefactor 
here results by expressing wave vectors in units of kp, 
energies in units of Ep, and radial positions in units of 
the Thomas-Fermi radius Rf = [2Ep j '(mw 2 )] 1 / 2 where 
ujq = (ujxUyUJz) 1 / 3 is the average trap frequency. 

Finally, to obtain an expression that can be directly 
compared with the experimental EDC spectra, it is suffi- 
cient to express the frequency Cj in Eq. (J3j) in terms of the 
single-particle energy E s = k 2 /(2m) — Cj via the relation 
£(fc;r) — uj = E 8 — /i(r), where = fi — mus^r 2 /2 is 
the local chemical potential in the harmonic trap. This 
yields eventually: 

48k 2 f°° 

EDC(fc, £,) = —-/ dr r 2 A(k, E s - fx(r))f(E s - M (r)). 
71 Jo 

(4) 

The numerical results obtained from Eq.Q are then con- 
voluted by a Gaussian broadening with a rms of about 
0.25Ep, corresponding to the experimental resolution. 

When the interparticle interaction is switched off, 
A(k,E s - n(r);r) is given by 5(E S - k 2 /(2m)). This 
defines the zero of the single-particle energy in the EDC 
curves as the energy of an isolated atom at rest. The 
chemical potential has thus disappeared from the free- 
particle branch k 2 /(2m), which remains positive for all k 
and can be used to reckon the value of the pscudogap. 

Experimental procedures 

We refer to Rcf. [5] for a detailed description of the ex- 
perimental techniques and procedures. Here, we add a 
few comments that are specifically relevant to the data 
presented in the main paper. 

These data arc taken at T/T c = 1.0 ± 0.1, where T c 
is determined in the trapped system by the vanishing of 
the measured condensate fraction. Note, however, that, 
because the density of the trapped gas is spatially inho- 
mogeneous, the local critical temperature decreases away 
from the cloud center. 

It was already remarked in the main paper that an ab- 
solute comparison can be made between the experimen- 
tal and theoretical EDC curves at given coupling and 
wave vector, in such a way that only their overall inte- 
gral over k and E s (and not the individual EDC curves) 
are normalized to unity. This is possible because the ex- 
perimental radio frequency data (from which the EDC 
curves are obtained) were taken for the first time over a 
wide range of w, such that their long high-frequency tail 
could be determined and accurately compared with the 
cj~ 3 / 2 behavior predicted theoretically 

The experimental data reported in Fig. 2(a) of the 
main paper have been analyzed in terms of the BCS-likc 
dispersion: 



E s (k) = /i - 
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FIG. 5. x 2 ~ fit to the experimental data. The experimental 
data of Fig. 2(a) of the main paper (circles) are fitted over 
the interval 0.0 < k/kF < 2.0 by the BCS-like dispersion of 
Eq.© (full lines). The free-particle dispersion k 2 /(2m) is also 
shown for comparison (thin full lines). 



where k^ ~ kp is the special wave vector about which 
the back-bending occurs and jl accounts for an overall 
(upward) displacement of the dispersion curves. Note 
that, contrary to the homogeneous case, in a trap fl is 
not related to the value of the thermodynamic chemical 
potential close to T c . 

A x 2 — analysis of the data in the interval 0.0 < k/kp < 
2.0 yields the fits shown in Fig. [5] The five values of kh 
thus obtained have been reported (together with the cor- 
responding error bars) in Fig. 3(b) of the main paper 
(there, the additional value for the coupling (kpap)^ 1 = 
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1.1 has been inferred from the experimental data re- 
ported in Ref.[H)- 

The same analysis also shows that pairs of (jl, A) with 
fi — A = constant produce comparable x 2 tests. From 
Eq.([5]) we note that p, — A = E s {k£) = E max corre- 
sponds to the maximum value of E s {k). For the five cou- 
plings here considered we obtain the values E max /Ep = 
(0.40, 0.24, -0.5, -1.1, -1.8), in the order. 

Using these values, one can extract a rough estimate 
of a (trap averaged) pseudogap energy, by relating them 
with the free-particle dispersion k 2 /(2m) at Ul- The k 2 - 
dispersion can, in fact, be considered as a lower bound 
to the dispersion of the upper branch in the EDC curves, 
which results from the two-peak structure of A(fc, w) in 
the presence of a pseudogap [1| and behaves like k 2 /(2m) 
for fc^ <C fc. [In the analysis of the experimental data 
the visibility of this upper branch is suppressed by the 
presence of the Fermi function.] The values we obtain 
for [kl/(2m)-E max ]/2 are (0.38, 0.34, 0.58, 0.82, 1M)E F 
for the five couplings of Fig. [5l which are in line with 
the expected trend for the pseudogap of a homogeneous 
system (cf. Fig. 17 of Ref.flj]). 

Determination of fci for a homogeneous system 

A comment is in order about the procedure for iden- 
tifying the Luttinger wave vector k^ for a homogeneous 
system, as reported in Fig. 3(a) of the main paper. 

Quite generally, A(k 7 uj) given by Eq.(p} has a pro- 
nounced peak when the following condition is satisfied 
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FIG. 6. Evolution in temperature of the dispersion oj(fc) ob- 
tained by following the peak at lower energy in A(k,uj) for a 
homogeneous system at unitarity. The four curves correspond 
to temperatures T/T c = (1.0, 1.2, 1.65, 4.0) from bottom to 
top. 
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£ k -ReS(k,w) = 



(6) 




and provided ImS(k, ui) 
like situation we write: 

S(k,w) 



is sufficiently small. In a BCS- 



A 2 

pg 



(7) 



with rj = + and Sfi = where fiL — k\/{2m). Note 

that the shift 6/j, is the part of the self-energy which is 
responsible for the persistence of a remnant Fermi surface 
about the (temperature dependent) Fermi wave vector of 
the underlying non-interacting system. A combination of 
Eqs.© and then yields: 
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pg 
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for the 
weight for k 



lower 

< 



2m 

(8) 

branch (which has the largest spectral 

kh), where oj is measured with respect to 
the chemical potential. The maximum value u w — A pg 
is for k = Ml where the back-bending occurs. For in- 
creasing temperature such that A pg closes up eventually, 
Eq.© yields accordingly k = k^ for lo = 0, which corre- 
sponds to the familiar condition for a Fermi liquid 0. 

Figure [5] shows a typical temperature evolution of the 
dispersion w(k) obtained by following the peak at lower 



FIG. 7. The dispersion (two arcs drawn by full lines) of 
the low-frequency peak of A(k, cu) for a homogeneous system 
with the density of the shell at r — r max for the coupling 
(^FtSF) 1 = 0.15, is compared with the corresponding dis- 
persion (dashed line) obtained multiplying A(k,u) by f(uj). 
A BCS-like fit to the two arcs is also shown (dotted line). 



energy in A(k,uj), from which the value of k = kh is 
extracted according to the above criterion. At sufficiently 
high temperatures when the dispersion w(k) crosses zero, 
kr, is seen to quickly converge to the value associated 
with the temperature dependent chemical potential of 
the non-interacting Fermi system. 

Additional theoretical analysis of the 
experimental data 

In the expression ((J) the presence of the Fermi func- 
tion f(cj) may be of considerable help for the analysis 
of the dispersion of the low-w peak, in situations when 
two broad non-Lorentzian structures in A(k,uj) merge 
together over a limited range of fc. This is because multi- 
plication of A(k, u) by f(ui) in that expression acts effec- 
tively as a "filter" for the low-w structures of A(k, uS), in 
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FIG. 8. The experimental EDC (circles) for the two cou- 
plings 0.45 and 0.57 are reproduced from Fig. 1 of the main 
paper, and compared with theoretical calculations (full lines) 
in which the temperature in the Fermi function has been de- 
creased to 0.7T C . The theoretical curves reported in Fig. 1 of 
the main paper are also reproduced here (dashed lines) . 



particular for those values of u> through which the back- 
bending occurs in the dispersion. 

This is shown explicitly in Fig. [71 where the dispersion 
of the low-w peak of A[k,uj) (corresponding to a homo- 
geneous system with the density of the shell at r max for 
the trap coupling (fcpOF) 1 = 0.15) is drawn (full line) 
only for those values of k for which two peaks in A(k, u>) 
appear clearly distinguishable. This procedure results in 
two arcs separated by an empty window. A single BCS- 
like fit (dotted line) to these two disconnected arcs via 
Eq.© provides the value A(r max )/Ep = 0.77, in reason- 
able agreement with the value determined for the whole 
trap. Figure [Jj shows also the dispersion (dashed line) 
obtained by multiplying A(k,u>) by f(cj), in such a way 
that the low-w peak can be smoothly followed even in 
the A;-window that had to be excluded before. This pro- 
cedure does not appreciably alter the values obtained by 
the BCS-like fit. 

This conclusion is consistent with the fact that the lack 
of a spectral depression in A(k, ui) in a limited range of k 
does not necessarily lead to disappearance of the pscudo- 
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T/T F 

FIG. 9. Specific heat per particle of a non-interacting Fermi 
gas (full line) reported over an extended temperature range. 
The dashed line corresponds to the linear behavior that holds 
when T/T F < 1. 



gap in integrated quantities, like the single-particle den- 
sity of states (DOS) [1], where a spectral depression sur- 
vives at much higher temperature than in A{k ks kp,u>) 
[cf. the inset of Fig. 4(b) of the main paper, where a cal- 
culation of the DOS is explicitly reported]. 

The presence of the factor f(u>) in Eq.@ obviously af- 
fects more the large-w than the low-w peak of the EDC 
curves. The discrepancies that are evident in the large- 
ly peak from Fig. 1 of the main paper, when comparing 
experimental and theoretical EDC curves at T c for the 
couplings 0.45 and 0.57, can accordingly be attributed to 
the larger absolute values of T c at which the theoretical 
spectra are calculated 0] , with respect to the experimen- 
tal values of T c at which the data are taken. 

In Fig. [8] we reproduce the experimental EDC (circles) 
from Fig. 1 of the main paper for the two couplings 0.45 
and 0.57, and compare them with the theoretical calcu- 
lations (full lines) in which the temperature in the Fermi 
function has been decreased to 0.7T C while the tempera- 
ture in A(k, lu) is kept at T c . This procedure is consistent 
with the fact that in this coupling range the theoreti- 
cal approach overestimates the absolute value of T c by 
about 30%, while close to T c the spectral function de- 
pends essentially on the relative temperature T/T c . This 
procedure, albeit empirical, goes in the right direction 
of reducing the height of the high-w peak of the EDC 
curves making it closer to the experimental value, while 
affecting only slightly the \ow-lu part of the EDC curves. 

The numerical difference between the theoretical and 
experimental values of T c for a homogeneous Fermi gas 
at unitarity is also evident from Fig. 4(c) of the main 
paper, although this difference is immaterial to the sake 
of the argument that was there raised. 

Specific heat of a non-interacting Fermi gas 

In Fig. 4(c) of the main paper the behavior of the spe- 
cific heat of a non-interacting Fermi gas was reported for 
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comparison over the temperature interval 0.2 < T/Tp < 
0.6, which was relevant to the experimental data shown 
in the same figure. 



It is interesting to show the same quantity over a more 
extended temperature range which reaches T = 0. This 
is done for < T/Tp < 1 in Fig. |H1 where the specific 
heat per particle of a non-interacting Fermi gas is re- 
ported vs T/Tp (full line) and compared with its linear 
approximation (£;£7r 2 /2) T/Tp (dashed line) that holds 
when T/Tp <C 1. Note that for T/Tp = 0.2 this linear 
approximation deviates from the full calculation already 
by about 20%. 



An analogous linear behavior is known to result for a 
Fermi liquid when T/Tp -C 1, although with a different 
slope reflecting the renormalization of the mass Q . 
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